J/ip suppression at y/s = 200 GeV 
in the comovers interaction model 



>n 

O ■ A. Capella a ), E. G. Ferreiro b ) 

o ■ 



in 



a) Laboratoire de Physique Theorique* 
Universite de Paris XI, Batiment 210, F-91405 Orsay Cedex, France 



b) Departamento de Fisica de Particulas, Universidade de Santiago de Compostela, 

<N ' 

■ 15782 Santiago de Compostela, Spain 

O 

in 
o 

>n ■ Abstract 
O 



Q H ' The yield of J/tp per binary nucleon-nucleon collision in AuAu and CuCu 

i- 



qj | collisions at = 200 GeV is computed in the framework of the dual parton 



model, supplemented with final state interaction (comovers interaction). For 
^ ' the latter we use the same value of the cross-section, a co = 0.65 mb, which 

describes the anomalous J/ip suppression observed at CERN-SPS energies. 
Several possibilities for the value of the absorptive cross-section are consid- 
ered. Shadowing is introduced in both the comovers and the J/ip yields. A 
comparison with the results at CERN-SPS, including a prediction for Inln 
collisions, is also presented. 

LPT-Orsay 05-23 
April 2005 



*Unite Mixte de Recherche UMR n° 8627 - CNRS 



1 



I. INTRODUCTION 



J I ip production in proton- nucleus collisions is suppressed with respect to the charac- 
teristic A 1 scaling of lepton pair production (Drell-Yan pairs). This suppression is generally 
interpreted as a result of the multiple scattering of a pre-resonance c — c with the nucleons 
of the nucleus (nuclear absorption). In these interactions, the c — c pair can transform into 
another c — c pair with vanishing projection into J/ip. The corresponding cross-section a a i JS 
is called absorptive cross-section. This interaction is generally described in the framework 
of a probabilistic Glauber model. However, at high energies, the coherence length increases 
and the projectile interacts with the nucleus as a whole. As a consequence the probabilistic 
Glauber formula breaks down [1] [2] and, thus, the extrapolation from CERN-SPS energies 
to RHIC ones (^Js = 200 GeV) is not straightforward. 

The NA50 collaboration has observed [3] the existence of anomalous J/ip suppression in 
Pb — Pb collisions, i.e. the J/t/j suppression in central Pb — Pb collisions clearly exceeds 
the one expected from nuclear absorption. Such a phenomenon was actually predicted 
by Matsui and Satz [4] as a consequence of deconfmement in a dense medium. It can 
also be described as a result of final state interaction of the c — c pair with the dense 
medium produced in the collision (comovers interaction). The final results [3] of the NA50 
collaboration can be described using an effective cross-section a co = 0.65 mb [5]. Since this 
is a low energy cross-section it is not expected to change significantly in going from CERN- 
SPS to RHIC energies. Therefore the prediction of the J/tp suppression due to comovers 
interaction at RHIC seems to be quite safe. However, at these energies it is necessary to 
introduce shadowing corrections. These are small at CERN-SPS energies. Moreover, they 
cancel to a large extend in the ratio of J/ip over Drell-Yan pair production measured by 
NA50. At RHIC energies, however, the Drell-Yan pair production is not measured and the 
J I ip suppression is presented as a ratio of the J / ip yield over the average number of binary 
nucleon-nucleon collisions - where the effect of shadowing is clearly present. 

II. THE MODEL 

The ratio R^(b) of the J/ip yield over the average number of binary nucleon-nucleon 
collisions, n(6), in AB collisions is given by 
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dN^(b)/dy 
n{b) 



dN^/dy 



jd 2 s a AB (b) n{b,s) S abs {b,s) S co {b,s) 
Jd 2 s cr AB (b) n(b,s) 



(1) 



Here a AB {b) = 1 — exp[— a pp ABT AB {b)} where T AB {b) = J d 2 sT A (s)T B (b — s) and T A (b) is 
the profile function obtained from Wood-Saxon nuclear densities [6] , and 



Upon integration of (2) over d 2 s we obtain the average number n{b) of binary nucleon-nucleon 
collisions at fixed impact parameter b. 

The factors S a b s and S co in (1) are the J/ip survival probability due to nuclear absorption 
and comovers interaction, respectively. 

In writing eq. (1) we have assumed that the J/ip yield in the absence of final state 
interaction (i.e. S abs = S co = 1) scales with the number of binary nucleon-nucleon collisions. 
In this case R AB coincides with the J / ip yield in pp collisions. 

a) Comovers interaction . The survival probability S co (b,s) of the J/ip due to co- 
movers interaction is obtained by solving the gain and loss differential equations which 
govern the final state interactions with the co-moving medium [7] 



where N J ^ and N co are the densities (i.e. number per unit of transverse surface) of J/ip 
and comovers (charged + neutral), respectively. In eq. (3) we have neglected a gain term 
resulting from the recombination of c and c into J/ip. This is natural in our approach since 
the cross-sections for recombination (gain) is expected to be substantially smaller than a co . 
The possibility of such a recombination, giving sizable effects at RHIC energies, has been 
considered by several authors [8] . It will be most interesting to see whether the data confirm 
or reject such an effect. 

In writing eq. (3) we have neglected transverse expansion and assumed a dilution 
in time of the densities due to longitudinal motion which leads to a r _1 dependence on 
proper time r. Eq. (3) can be solved analytically. The solution is invariant under 
the change r — > cr. Thus, the result depends only on the ratio T//Vo of final over ini- 
tial time. Using the inverse proportionality between proper time and densities, we put 
Tf/r = N co (b, s,y)/N pp (y), i.e. we assume that the interaction stops when the densities 
have diluted, reaching the value of the pp density at the same energy. At ^Js = 200 GeV 



n(b, s) = AB a pp T A (s) T B (b - s)/a AB (b) . 



(2) 



dN J ^(b,s,y) 



a co N J ^(b,s,y) N co (b,s,y) 



(3) 



dr 
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and y* ~ 0, N pp (0) = \(dN ch /dy)^ =0 /nR 2 p ~ 2.24 fin" 2 . Note the increase in the pp density 
form 1.15 fm~ 2 at CERN-SPS to 2.24 fm -2 . Since the corresponding increase in the AuAu 
densities is approximately the same, the average value of Tf/r is about the same at the two 
energies (of the order of 5 4- 7) . 

The solution of eq. (3) is given by [5] 

S co {b, s) = N J/Winal) {b, s, y) /N J/ ^ initial) {b, s, y) 

= exp [-a co N co (b, s, y)£n (N co (b, s, y)/N pp (0))} . (4) 

b) Comovers density in the dual parton model . The main ingredient in order to 
compute the survival probability S 00 is the comovers density N co . Note that N co is the 
comovers density at initial time r , i.e. the density produced in the primary collision. It can 
be computed in the dual parton model [9]. It can be seen from eq. (6.1) of ref. [9] that this 
density is given by a linear combination of the average number of participants and of binary 
nucleon-nucleon collisions, and for A = B can be written as 

3 HN ch 3 

Nf? s (b,8,y) = - -^(b,s,y) = - [d(&) n A (b, s) + C 2 (b) n(b,s)\ (5) 

where 

n A (b, s) = A T A (s) [1 - exp (-a pp B T B (b - s))}/ a AB (b) (6) 

and n(b,s) is given by eq. (2). The factor 3/2 takes care of the neutrals. The coefficients 
Ci(b) and C^fr) are obtained from string multiplicities which are computed in DPM as 
a convolution of momentum distributions functions and fragmentation functions. These 
functions are universal, i.e. the same for all hadronic and nuclear processes. Thus, we use 
the same expressions as at CERN energies. For details see [10]. The numerical values of 
Clip) and C 2 (6) in AuAu and CuCu collisions at y/s = 200 computed in the rapidity interval 
\y*\ < 0.35 for various values of b and per unit rapidity are given in Table 1. We also give 
in this table the corresponding values for PbPb and Inln at pi ab = 158 GeV/c. 

We see from Table 1 that C 2 is significantly larger than C\ at RHIC energies. Thus, 
DPM leads to multiplicities which have a behaviour closer to a scaling with the number of 
binary collisions rather than to a scaling with the number of participants. Actually, with 
increasing energies the ratio CijC\ increases and one obtains a scaling in the number of bi- 
nary collisions. This is a general property of Gribov's Reggeon Field Theory which is known 
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as AGK cancellation [11] - analogous to the factorization theorem in perturbative QCD and 
valid for soft collisions in the absence of triple Pomeron diagrams. It is well known that 
this behaviour is inconsistent with data which show a much smaller increase with centrality. 
As discussed in detail in [10] such a discrepancy is due to shadowing which is important at 
RHIC energies and has not been taken into account in eq. (5). This is precisely the meaning 
of label NS (no shadowing) in this equation. 

c) Shadowing . Following ref. [10,12,13] shadowing corrections are computed, without 
free parameters, in terms of the measured value of the diffractive cross-section. Indeed, in 
the framework of Gribov's Reggeon Field Theory, the same triple Pomeron diagrams which 
describe high mass diffraction are responsible for the shadowing corrections. While the con- 
tribution of the triple Pomeron diagram to high mass diffraction is positive, its contribution 
to the total cross-section is negative, due to the presence of s-channel discontinuities which 
correspond to interference terms. Thus, the triple Pomeron diagrams (with triple Pomeron 
coupling determined from high-mass diffraction data) produce a decrease of the charged 
yield as given by eqs. (5)-(7), thereby violating the AGK cancellation. In AB collisions this 
reduction is given''' by [10,12] 

1 1 

^ *' V) = I + A F h (y) T A (s) 1 + B F h (y) T B (b - s) ' (?) 

Here the function F(y) is given by the integral of the ratio of the triple Pomeron cross-section 
d 2 a ppp jdYdt at t — over the single Pomeron exchange cross-section a P : 

rY max 1 d 2 a PPP 



f J max 

F h (y) = 4tt / dY 

J T'rn in 



op dYdt 



= C [exp (AY max ) - exp (AY min )} (8) 

t=o 

where Y min = £n(RAm N / \f?>) , A = 0.13 and C = 0.31 fm 2 . The value of Y max depends on 
the rapidity of the produced particle. For y* = we have Y max = l/2£n(s/rri^) where mr is 
the transverse mass of the produced particle. For charged particles we use uit = 0.4 GeV 
and for a J/ip m T = 3.1 GeV. 



1"A more sophisticated calculation using other triple Regge diagrams, with parameters constrained 
from HERA data can be found in ref. [13]. The results, however, are very similar to the ones 
obtained from eqs. (7)-(8). 
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It has been shown in [10] that with the shadowing resulting from eqs. (7)- (8) a good 
description of the centrality dependence of charged multiplicities at mid-rapidities is obtained 
at RHIC energies (y/s = 130 and 200 GeV). More precisely one has 

N c °(b, s, y) = N c N ° s (b, s, y) S*(b, s, y) (9) 

where the two factors in the r.h.s. are given by eqs. (5) and (7), respectively. 

With this expression of the density of comovers we can compute the J/ip survival prob- 
ability S co , eq. (4). The J/ip suppression R J /^ is given by eq. (1) with the following 
replacement in its numerator 

n(b,s)^n(b,s) S^(b,s,y) . (10) 

Indeed, as discussed above, in writing the numerator of eq. (1) we have assumed that the 
J/ijj yield in the absence of final state interaction (S abs = S co = 0) scales with the number of 
binary collisions. This is only true when shadowing is neglected. The effects of shadowing 
on the J/ip yield are introduced with the replacement (10) in the numerator of eq. (1). 



d) Nuclear absorption . The formula for nuclear absorption used in the literature is 
obtained in a probabilistic Glauber model. One has 

qab s (h \ _ [1 - exp(-A T A (s) a abs )}[l - exp(-ff T B (b - s)a abs )] 
S ABT A (s)T B (b-s) • 

As discussed in the Introduction, this formula breaks down at high-energy due to the increase 
of the coherence length [1] [2]. In the limit of s — > oo, the relevant equation is quite simple. 
The change consists in the replacement 

1 



'l/a ahs ) [1 - exp (-Gabs A T A (b))\ =>- A T A (b) exp 



2 : ct C c_at A T A (b) 



(12) 



in each of the two factors in the numerator of (11). The corresponding formula at finite 
energies which interpolates between (1) and (2) has also been derived [1]. The change in 
going from (11) to (12) is twofold. There is a change in the form of the expression and, 
moreover, a abs has been replaced by the total cc — N cross-section <j C c-n- If &cc-n ~ &abs the 
change from low energies to asymptotic ones is small. Indeed the two expressions coincide 
at the first and second order in the development of the exponential, and since a a b s is small, 
the low energy result will not be significantly changed. However, if <7 C c-n is significantly 
larger than its absorptive part, cr abs , the J/ip suppression due to final state interaction 
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within the nucleus will be larger at high energies. The latter possibility seems to be ruled 
out by preliminary data [14] on dAu collisions which show a rather small suppression at 
mid-rapidities. 

In the next section we present the calculation of J/tp suppression in AuAu and CuCu 
collisions using eq. (11) with the value a a bs = 4.5 obtained from the pA data at CERN-SPS. 
With these values of a a bs the results obtained with eqs. (1) and (2) are practically the same. 
We also present the results obtained with smaller values of cr a b s {(Tabs = 0, 1 mb and 3 mb). 



At mid-rapidities, this quantity can be computed at various centralities using the coefficients 
Ciib) and C2(b) in Table 1 and eq. (7). The calculations for AuAu collisions at mid-rapidities 
are shown in Fig. 1 and compared with experiment [15]. As predicted in [10] a reasonable 
description of the data is obtained. An increase by a factor 1.13 between y/s = 130 GeV and 
y/s = 200 GeV for central collision was also predicted in [10] - in agreement with present 
data. 

In Fig. 2 we compare the result of our calculations of Rab(^) i n ec f (1) f° r different 
systems : PbPb* and Inln at CERN-SPS (pi a b — 158 GeV/c) and AuAu and CuCu at 
yfs = 200 GeV. In all cases the normalization is arbitrary but the same for all. It corresponds 
to taking dNij^ /dy = 1 in eq. (1). Also in all cases we use a co = 0.65 mb and a abs = 4.5 mb. 
An important feature of our results is that, at a given energy, the results for the lighter 
systems are rather close to the ones for the heavier ones, at the same values of N part . We 

^The results for PbPb are identical to those in the first paper of [5], except that in [5] the ratio 
J /if) over DY was plotted versus Et (the energy deposited in the NA50 calorimeter). Moreover, for 
large Et (beyond the "knee" of the Et distribution) the effect of the fluctuation in the comovers 
multiplicity was included. This is not the case in Fig. 2 since, in a plot versus N part , such a situation 
does not arise. 



III. NUMERICAL RESULTS 



We compute first the inclusive charged particle multiplicity given by 




(13) 
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also see that the J/ip suppression is much larger at RHIC energies and reaches a factor 10 
for central AuAu collisions. 

In Fig. 3 we present again the result of our calculation of RauaJP) m ec l- w ^ n the 
normalization given by the measured value of dN^/dy at = 200 GeV [16]. All curves 
are obtained with a co = 0.65 mb and different values of a abs (cr abs = 0, lmb, 3 mb and 4.5 
mb). The suppression for central collisions varies between a factor of 6 for a abs = and a 
factor of 10 for a abs = 4.5 mb. Even in the former case the suppression is twice as large as 
the one obtained in a QCD based nuclear absorption model [17]. 

The results in Fig. 2 for PbPb and Inln has been obtained without including shadowing. 
Using eq. (7) it turns out that, at CERN-SPS, the shadowing on the J/ip is negligibly small. 
However, for the comovers (m T = 0.4 GeV), its effect is of the order of 15 %. Here the values 
of Y max and Y min in eq. (8) are quite close to each other and our equations (and in particular 
the expression of Y min ) are not accurate enough for a reliable calculation. If, however, an 
effect of shadowing of the order of 15 % is present, the values of C\ and C 2 in Table 1 should 
be increased by the same amount in order to restore agreement [10] with the experimental 
values of the charged multiplicities in PbPb at pi ab = 158 GeV/c. This, in turn, would result 
in a larger J ftp suppression at RHIC. The maximal effect occurs in the case a abs = and is 
of the order of 20 % for central AuAu collisions. 

IV. CONCLUSIONS 

In a comovers interaction framework we have computed the yield of J/ip per binary 
nucleon nucleon collision versus the number of participants in PbPb and Inln collisions 
at CERN-SPS (p lab = 158 GeV/c) and in AuAu and CuCu at ^ = 200 GeV. At RHIC 
energies shadowing corrections to both the J/ip and the comovers multiplicities are very 
important and have been included in the calculations. We have found that, at a given 
energy, the J/ip suppression for the lighter and heavier systems are similar, at the same 
value of N part . We have also found that the J ftp suppression at RHIC is significantly larger 
than at SPS. For central AuAu collisions it reaches a factor of 10 for a abs = 4.5 mb and a 
factor 6 for a abs = 0. The value of a abs has to be determined from the dAu data. Preliminary 
results [16] favor a rather small value, a ab m 1 mb. 

We have argued that these values could be underestimated by about 20 %. Experimental 
values of the J/ip suppression significantly smaller that the one in Fig. 3 would not be 
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consistent with the comovers interaction model, at least in its present formulation. 

Finally, an important difference between the J/tp suppression pattern in a comovers 
interaction model and in a deconfining scenario is that, in the former case, the anomalous 
supression sets in smoothly from peripheral to central collisions - rather than in a sudden 
way when the deconfining threshold is reached. The NA50 results have not allowed to 
disentangle these two possibilities. However, at RHIC energies, the relative contribution of 
the comovers is strongly enhanced in our approach, and a clear cut answer to this important 
issue should be obtained. 
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TABLES 



b 


Ci AuAu 


C 2 AuAu 


d CuCu 


C 2 CuCu 


d PbPb 


C 2 PbPb 


Ci Inln 


C 2 Inln 


0. 


1.0274 


1.7183 


1.0330 


1.8196 


0.7102 


0.3975 


0.7480 


0.4312 


1. 


1.0276 


1.7206 


1.0334 


1.8239 


0.7115 


0.3987 


0.7485 


0.4317 


2. 


1.0278 


1.7228 


1.0338 


1.8320 


0.7152 


0.4020 


0.7527 


0.4357 


3. 


1.0286 


1.7340 


1.0342 


1.8437 


0.7208 


0.4070 


0.7599 


0.4428 


4. 


1.0293 


1.7448 


1.0347 


1.8592 


0.7283 


0.4136 


0.7696 


0.4526 


5. 


1.0302 


1.7574 


1.0352 


1.8787 


0.7376 


0.4218 


0.7810 


0.4646 


6. 


1.0310 


1.7722 


1.0357 


1.9014 


0.7488 


0.4320 


0.7945 


0.4793 


7. 


1.0320 


1.7908 


1.0361 


1.9258 


0.7617 


0.4445 


0.8112 


0.4985 


8. 


1.0330 


1.8121 


1.0364 


1.9505 


0.7764 


0.4597 


0.8290 


0.5198 


9. 


1.0340 


1.8374 


1.0364 


1.9754 


0.7929 


0.4776 


0.8475 


0.5430 


10. 


1.0349 


1.8665 


1.0363 


2.0006 


0.8112 


0.4985 


0.8664 


0.5681 


11. 


1.0357 


1.8990 


1.0360 


2.0259 


0.8308 


0.5220 


0.8855 


0.5949 


12. 


1.0362 


1.9308 


1.0356 


2.0515 


0.8503 


0.5466 


0.9046 


0.6235 


13. 


1.0364 


1.9580 


1.0349 


2.0772 


0.8673 


0.5698 


0.9233 


0.6536 



TABLE I. Values of C\ and C 2 in eq. (5) as a function of the impact parameter b. The second 
and third columns correspond to AuAu collisions and the forth and fifth to CuCu collisions both 
at yfs = 200 GeV. The values, calculated in the range —0.35 < y* < 0.35, are given per unit 
rapidity. The following columns refer to PbPb and Inln at pi a b = 158 GeV/c and are computed 
in the rapidity range of the NA50 dimuon trigger < y* < 1. 
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FIGURES 
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FIG. 1. Multiplicity of charged particles per participant pairs, versus N part , computed from 
(13) are compared to experimental data from PHENIX [15]. 
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FIG. 2. R J Js(b) for AuAu collisions at ^/s = 200 GeV (full curve), CuCu collisions at yfa = 200 
GeV (dashed curve), PbPb at pi a b = 158 GeV/c (dotted curve) and Inln at pi^ = 158 GeV/c 
(dashed-dotted curve). All the results have been obtained with a co = 0.65 mb and a a bs = 4.5 mb. 
The normalization, the same for all four curves, is arbitrary. It corresponds to taking dNpp^ jdy = 1 
in eq. (1). 
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FIG. 3. R J j[fi{b) for AuAu collisions at \fs = 200 GeV multiplied by the dilepton branching 
ratio, normalized to the measured value in pp collisions [16]. From up to down: result with 
a co = 0.65 mb and a a b s = mb (dashed curve), result with a co = 0.65 mb and a a bs = 1 rnb 
(dotted-dashed curve), result with o co = 0.65 mb and a a bs = 3 mb (dotted curve) and result with 
a co = 0.65 mb and a a b s = 4.5 mb (full curve). 
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